Abstract. Using the sample of radio selected BL Lacer- 
tae objects (RBLs) and X-ray selected BL Lacertae ob- 
jects (XBLs) presented by Sambruna et al. (1996), we 
calculated the luminosities of radio, optical and X-ray of 
each source and made the statistical analysis among the 
luminosities at different wave-bands, broad-band spectral 
indices from radio to X-ray (a rx ) and peak frequencies 
{v p ). Our results are as follows: (i) there is a positive cor- 
relation between radio luminosity L r and a rx and a nega- 
£N| ' tive correlation between L r and v v . High-energy peak BL 
. Lacs (HBLs) and low-energy peak BL Lacs (LBLs) can 
be distinguished very well, the dividing lines are probably 
those of logi r = 43.25 (erg/sec) and a rx >(or <)0.75 for 
L T - a rx plot and those of logi r < 43.25 (erg/sec) and 
logfp > 14.7 for the L r - v p plot; (ii) there is a weak pos- 
itive correlation between optical luminosity L Q and a rx 
and a negatively weak correlation between L a and is p ; (iii) 
there is no correlation between X-ray luminosity Lx and 
a rx or between Lx and v v . From our analysis, we find that 
synchrotron radiation is the main X-ray radiation mecha- 
nism for HBLs while inverse Compton scattering for LBLs. 
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Most BL Lac objects have been identified in either 
radio survey ( radio-selected BL Lac objects (RBLs) or 
X-ray survey (X-ray-selected BL Lac objects (XBLs)). 
These two kinds of BL Lac objec ts exist many signifi- 
cant differences( for example see, Ledden, O'Dcll, 1985 ; 
Ghiscllini et a!.1986t [Stocke et al. 199l| ). The main differ 



ences between RBLs and XBLs are as follows. First, there 
is a difference in spectral energy distributions (SEDs), 
RBL-like objects show a IX > 0.75, while XBL-like 
objects show q rx < 0.75 (see, e.g., Ledden, O'Dell, 1985; 



Stocke et al. 1985 



Stocke et al. 1991 



Giommi et al. 199C 



Schachter et al. 1993; Giommi & 



Padovani 1994). It should be pointed out that Padovani 
& Giommi (1995) introduced the distinction between 
high-energy peak BL Lacs (HBLs) and low-energy peak 
BL Lacs (LBLs), for objects which emit most of their 
synchrotron power at high (UV-soft-X) and low (far-IR, 
optical) frequencies respectively. However, a quantitative 
distinction can be drawn on the basis of the ratio between 
radio and X-ray fluxes. It can also use the broad-band 
spectral index a lx to distinguish HBLs and LBLs: 
a rx < 0.75 for HBLs and a rx > 0.75 for LBLs (Giommi 
& Padovani 1995). Secondly, there is a difference in 
the peak frequency dis tributi ons ( Giommi et al. 1995 ; 
Padovani fc Giommi, 1995 ; Urry fc Padovani, 1995 ; 

Padovani fc Giommi, 1996 ; Lamer et al. 1996 ). Sam- 
bruna et al.(1996) made a parabolic fit to multi-band 
flux data to obtain the peak frequency of the synchrotron 
emission, v Vl of sources, when studying quantitatively 
how the peak-emission frequency of the synchrotron 
emission can be used to distinguish RBL-like and 
XBL-like objects. Making use of the peak frequency, 
Up, of the synchrotron emission estimated by Sambruna 
et al.(1996). Qin et al.(1999) found that, in the four 
different regions divided by the arx = 0.75 line and the 
logz/p = 14.7 line, all RBL-like objects lie in the upper 
left region, while most XBL-like objects are within the 
lower right region, with only a few sources being located 
in the lower left region and no sources being located in 
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the upper right region. Considering the proper physical 
origin of the X-ray emission for different classes of BL 
Lac objects, we calculated the peak frequency v v of 
the synchrotron emission of each source in sample of 
Sambruna et al. (1996). It was found that all RBL-like 
objects, defined by a rx > 0.75, are located in the 
upper-left region, while all XBL-like objects, defined by 
Q^rx < 0.75, are inside the lower-right region. No sources 
are in the lower left and upper-right regions, suggesting 
that the classificatory criteria in terms of the peak 
frequency should be \o%v v — 14.7 and the classificatory 
criteria of the a rx is equal to the classificatory criteria 
of the \ogv p (Dong, Mei, & Liang, 2002). This provides 
evidence supporting what Giommi et al. (1995) proposed: 
RBL-like and XBL-like objects can be distinguished by 
the difference in the peak frequency of the synchrotron 
emission. 

Sambruna et al.(1996) found that the position of the 
peak frequency log v p is linked to the luminosity. For more 
luminous objects, the peak of the synchrotron power lo- 
cate at lower frequencies. Fossati et al.(1998) studied the 
spectral energy distributions of three subclass samples 
of blazars. They also found that despite the differences 
in continuum shapes for sub-classes of blazars, a unified 
scheme is possible, blazar continua can be described by 
luminosity as a fundamental parameter. 

In this paper, we will study how the luminosity of each 
energy band can be used to distinguish RBLs and XBLs. 
We calculated luminosity of each energy band for RBLs 
and XBLs presented in Table 1 of Sambruna et al.(1996) 
in section 2. We present our analysis results in section 3, 
and give our discussion and conclusions in section 4. 



1. Luminosity of BL Lac objects 

We calculated the luminosities for the sample of RBLs 
and XBLs presented by Sambruna et al.(1996). We list 
the names, red-shifts, broad-band spectral indices, the en- 
ergy fluxes at radio, optical and X-ray bands and peak 
frequencies of 23 XBLs in Table 1 and 29 RBLs in Table 

2, respectively. Using the logrithmic parabolic form given 
by Laudau et al (1986) and the data of the energy fluxes 
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at radio, optical and X-ray bands, Dong et al. (2002) cal- 
culated the peak frequencies of 23 XBLs in Table 1 which 
are used here. Sambruna et al.(1996) made a parabolic fit 
to the radio, millimeter, IR, optical, and X-ray fluxes data 
to obtain the peak frequency, v p , of the spectral energy, 
where, for some XBLs, X-Ray data were not included in 
the fit. The data of peak frequency v pi of the synchrotron 
emission in Table 2 are made use of Sambruna et al.(1996) 
except for objects 1652+398 and 2005-489, because the 
both sources are HBLs (a rx < 0.75). We use the v p data 
of two sources calculated by Dong, Mei & Liang (2002). 
Therefore, in our sample, there are 27 LBLs which con- 
sists of RBLs except for 1652+398 and 2005-489 in Table 
2 and 25 HBLs which consists of 23 XBLs and two RBLs 
(1652+398 and 2005-489). 

We calculated rest-frame luminosities of radio, optical 
and X-ray energy band using the fluxes of BL Lac ob- 
jects in Table 1 and Table 2. For the K-correction, the 
flux densities were multiplied by (1 + z) a , where a is 
the power-law spectral index in the appropriate energy 
band (F v oc v~ a ). For RBLs , we used a r — 0.2, and 
a op t = 1.05 ( Falomo et al. 1994), while for XBLs we used 
a r = and a op t = 0.65 (Falomo et al. 1994). We used indi- 



vidual spectral indices in X-ray provided by Sambruna et 
al.(1996). The redshifts of 0048-097,1147+245 and 1519- 
273 presented in table 2 of Sambruna et al.(1996) are not 
available. A mean value of z = 0.56 is accordingly adopted 
for these sources. The flux densities were converted to lu- 
minosities using Ho — 75 Km S" 1 Mpc" 1 and qo = 0.5, 
assuming isotropic emission. We calculated the radio lu- 
minosity L Xl optical luminosity L Q and X-ray luminosity 
L x in rest frame. 

2. Results 

We now consider the relations among radio luminosity, 
broad-band spectral index and peak frequency. We plot L r 
versus a rx in figure la, which shows that the radio lumi- 
nosity increases as broad-band spectrum index increases. 
A linear regression analysis yields 

logL r = (6.64 ± 0.49)a rx + (38.24 ± 0.37), (1) 
where the correlation coefficient is r = 0.886 and the 
chance probability is p — 2.26 x 10 -17 , n=52. Obviously, 
the correlation is very strong. From this figure, four dif- 
ferent regions are divided by the lines a rx = 0.75 and 
logi r = 43.25 (ergs -1 ), all LBLs are located in the upper- 
right region; all HBLs are inside the lower-left region(two 
RBLs, 1652+398 and 2005-489 are HBLs); no sources are 
in the lower-right region and the upper-left region. We plot 
L T versus v p in figure lb, which show that the radio lu- 
minosity L r decreases as the peak frequency v p increases. 
The correlation analysis gives 

logX r = -(0.67 ±0.06) log v p + (53.13 + 0.92), (2) 
where the correlation coefficient is r = —0.841 and the 
chance probability is p — 4.42 x 10 -14 , n=52. The cor- 
relation is also very strong. In the four different regions 



divided by the lines v p = 14.7 and log L r = 43.25, all LBLs 
are located in the upper-left region; all HBLs are inside 
the lower-right region(two RBLs, 1652+398 and 2005-489 
are HBLs) ; no sources are in the lower-left region and the 
upper-right region. 

In figure 2, we consider the correlation between the 
optical luminosity L a and a rx and that between L Q and 



Up. The results are 



(3) 



logL = (2.45 ± 0.45)a rx + (43.94 ± 0.33) 
with r = 0.614 and p = 2.70 x 10 -6 , n=52 and 
logL = -(0.28 ± 0.05) logvp + (49.93 ± 0.68) (4) 
with r = —0.660 and p = 2.52 x 10~ 7 , n=52 respectively. It 
can be seen that (i)there are a positive correlation between 
I/ and a rx and a negative correlation between L a and v p , 
and (ii)HBLs and LBLs can be distinguished by the line 
of a rx w 0.75 (see Fig. 2a) or the line of v p » 14.7 (see 
Fig. 2b), but they can not be distinguished by optical 
luminosity. 

Finally, we study the relation of the X-ray luminosity 
I x with a rx and v p . The correlation analysis reads 
logi x = -(1.41 ± 0.43)a rx + (46.08 ± 0.32) (5) 
for L x versus a rx , where r = —0.42 and p = 0.002. 
\ogL x = (0.11 ± 0.05) \ogVp + (43.43 ± 0.73) (6) 
for L x versus v p , where r = 0.30 and p = 0.03. Obviously, 
there is no any correlation between L x and a rx as well 
as between L x and v p . We plot the results in Fig. 3. The 
distributions of HBLs and LBLs in plot of L x versus a rx ( 
Fig. 3a) and L x versus v p (Fig. 3b) are confusion although 
the lines of a x « 0.75 and v p ss 14.7 can be drawn. 

The results of the regression analysis are listed in Table 

3. 



3. Discussions and conclusions 

Sambruna et al.(1996) have made the statistical analy- 
sis of the relations of bolometric luminosity Lb (as esti- 
mated from parabolic fits to the multi-wavelength from 
radio to X-ray bands) with position of the peak frequency 
log u p and with the broad-band spectral index a TX for the 
blazars. They found that (i)there is a correlation between 
Lb and a rx at > 99% confidence, which may indicate that 
more luminous objects have steeper a rx , and (ii) more lu- 
minous sources have smaller peak frequencies, since the 
average bolometric luminosity increases from XBLs to 
RBLs to flat spectrum radio quasars (FSRQs), on aver- 
age, FSRQs are more luminous and have lower peak fre- 
quencies while XBLs are less luminous and have higher 
peak frequencies. In their analysis, the XBLs and RBLs 
are not distinguished very well in cither LB-a IX plot or 
Lb-u p plot. In this paper, we have studied the relations 
of radio, optical and X-ray luminosities with broad-band 
spectral index (a rx ) and peak frequency (u p ) for 52 BL Lac 
objects (including 23 XBLs and 29 RBLs). It is different 
from the analysis of Sambruna et al. (1996), we analyzed 
the relation of luminosity at the each energy band with 
a rx and v p , respectively. 
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From our analysis, we found that radio luminosity L T 
can be used to distinguish HBLs and LBLs. Our results 
indicate that there is a very strong positive correlation 
between L T and a rx and a negative correlation between 
L T and v p , more importantly, HBLs and LBLs can be dis- 
tinguished very well in both L Y -a IX plot (see Fig. la) and 
L T -i/ p plot (see Fig. lb). We have also analyzed the rela- 
tions of the optical and X-ray luminosities with cn rx and 
v p , respectively. We found that the optical luminosity has 
a good correlation with both a rx and v p (see Fig. 2a and 
2b), but X-ray luminosity shows no any correlation with 
a IX or v p (see Fig. 3a and 3b). Obviously, either L Q or L x 
cannot be used to distinguish HLBs and LBLs. Therefore, 
a rx , v p and L r are equivalent in the classifications of HBLs 
and LBLs: HBLs fall in the region divided by the lines of 
a rx < 0.75 (or log ^ p > 14.7 Hz) and log L T < 43.25 erg s _1 
while LBLs in the region divided by the lines of a rx > 0.75 
(or log v p < 14.7 Hz) and logi r > 43.25 erg s _1 . 

It has been known that the observed spectral energy 
distributions of BL Lac objects show a peak between IR 
and X-rays, and the possible radiation mechanism is syn- 
chrotron radiation. In Fig. la of L r and a rx plot, HBLs 
distribute in the region divided by the lines of a rx < 0.75 
and logi r < 43.25 erg s _1 . Using the definition of a rx 
and a rx < 0.75, we have logL x > logi r + 1 erg s _1 for 
HBLs where v r = 5 GHz and v x = 1 keV are used. Since 
logL r < 43.25 erg s _1 , we have logL x > 44.25 erg s _1 for 
HBLs while logL x < 44.25 erg s _1 for LBLs. Compared 
to the result shown in Fig. 3a, the observed X-ray lumi- 
nosities of HBLs in our sample are above the lower limit 
of log L x > 44.25 erg s _1 , indicating that X-rays for HBLs 
are produced by the synchrotron radiation. From Fig. 3a, 
the observed X-ray luminosities for most LBLs do not sat- 
isfy the upper limit of L x < 44.25 erg s _1 . We believe that 
inverse Compton scattering have more important contri- 
butions to X-ray emission from these LBLs. 

Finally, we would like to point out that we have made 
the analysis about the relation of energy fluxes at radio, 
optical and X-ray bands with a rx and v p . We find that 
radio energy flux can be distinguish HBLs and LBLs very 
well. 
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Table 1. Sample of X-ray selected BL lacertae objects (XBLs). 
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Name 


z 




a 


* 

rx 


(5 GHz) 


f?(mJy) 
(5500 A) 


(1 keV) 


(Hz) 


0112.1+0903 





.339 





,58 


0.0014 


0.047 + 0.001 


0.05 


15.62 


0158.5+0019 





.299 





,51 


0.0113 


0.21 + 0.06 


1.2 


17.34 


0205.7+3509 





.318 





,46 


0.0036 


0.10 + 0.005 


0.90 


18.06 


0257.9+3429 





.247 





,65 


0.010 


0.25 + 0.02 


0.10 


15.03 


0317.0+1834 





.190 





,58 


0.017 


0.36 + 0.09 


0.54 


15.79 


0419.3+1943 





.512 





,53 


0.008 


0.09 


0.75 


18.31 


0607.9+7108 





.267 


0. 


.70 


0.0182 


0.09 


0.07 


15.23 


0737.9+7441 





.315 





,56 


0.024 


0.64 


1.30 


16.09 


0922.9+7459 





.638 





,55 


0.0033 


0.044 ± 0.002 


0.21 


17.21 


0950.9+4929 





.207 





,51 


0.0033 


0.122 ± 0.04 


0.27 


16.13 


1019.0+5139 





.141 





,45 


0.0024 


0.22 


0.93 


16.69 


1207.9+3945 





.615 





,52 


0.0058 


0.10 


0.55 


17.38 


1221.8+2452 





.218 





67 


0.0264 


0.42 ± 0.09 


0.18 


15.00 


1229.2+6430 





.164 





,56 


0.042 


0.55 + 0.17 


2.05 


16.71 


1235.4+6315 





.297 





,55 


0.007 


0.14 + 0.02 


0.31 


16.16 


1402.3+0416 





.200 





,57 


0.0208 


0.88 + 0.37 


0.68 


15.39 


1407.9+5954 





.495 





,66 


0.0165 


0.07 + 0.01 


0.10 


15.76 


1443.5+6349 





.299 


0. 


.58 


0.0116 


0.06 


0.35 


17.71 


1458.8+2249 





.235 





,58 


0.0298 


1.01 + 0.20 


0.78 


15.39 


1534.8+0148 





.312 





.61 


0.034 


0.15 + 0.05 


0.74 


17.40 


1552.1+2020 





.222 





,54 


0.0375 


0.44 + 0.08 


2057 


17.45 


1757.7+7034 


0. 


.407 





50 


0.0072 


0.18 


0.92 


17.15 


2143.3+0704 





.237 





,61 


0.050 


0.32 + 0.04 


0.78 


16.24 



'These data are taken from Tables 1 and 2 of Sambruna et al. (1996). 
*These data are taken from Tables 1 of Dong,Mei& Liang (2002). 
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Table 2. Sample of radio selected BL lacertae objects (RBLs). 



Name 


z* 


a 


rx 


^r*(Jy) 


fo'(mJy) 




logzV 










(5 GHz) 


(5500 A) 


(1 keV) 


(Hz) 


0048-097 







.75 


1.110 ± 0.583 


2.41 ± 1.63 


0.77 


13.84 


0118-272 


>0.557 





.86 


1.145 ± 0.075 


1.92 ± 0.38 


0.20 


14.49 


0235+164 


0.940 


o 


.76 


1.81 ± 0.54 


1.44 ± 1.06 


1.56 


13.39 


0426-380 


> 1.030 





.90 


1.15 ± 0.03 


0.11 


0.09 


13.22 


0454+844 


0.112 


1. 


.01 


1.41 ± 0.14 


0.70 ± 0.36 


0.02 


13.81 


0537-441 


0.896 





.82 


3.93 ± 0.17 


1.49 ± 0.43 


0.78 


14.07 


0716+714 


>0.300 


o 


.75 


0.86 ± 0.18 


2.96 


1.17 


13.79 


0735+178 


>0.424 


o 


.88 


2.13 ± 0.50 


3.22 ± 1.56 


0.22 


14.03 


0814+425 


0.258 





.98 


1.86 ± 0.68 


0.26 ± 0.04 


0.05 


13.34 


0851+202 


0.306 


o 


.84 


2.99 ± 0.56 


6.08 ± 5.91 


0.70 ± 0.25 


13.72 


0954+658 


0.367 


o 


.88 


0.90 ± 0.38 


0.86 ± 0.24 


0.16 


14.09 


1144-379 


1.048 





.82 


1.61 ± 0.96 


0.62 ± 0.37 


0.41 


13.75 


1147+245 




o 


.92 


0.82 ± 0.12 


1.53 ± 0.36 


0.05 


14.58 


1308+326 


0.997 





.91 


2.26 ± 0.40 


2.23 ± 1.53 


0.13 


13.83 


1418+546 


0.152 





.85 


1.22 ± 0.38 


2.72 ± 0.82 


0.30 


13.85 


1519-273 







.86 


2.17 ± 0.25 


0.47 ± 0.35 


0.39 


13.17 


1538+149 


0.605 


o 


.93 


1.53 ± 0.42 


0.32 ± 0.10 


0.09 


13.56 


1652+398 


0.033 





.67 


1.27 ± 0.10 


15.65 ± 4.52 


8.30 


15.01 


1749+096 


0.320 





.92 


1.44 ± 0.36 


1.18 ± 0.54 


0.14 + 0.01 


13.27 


1749+701 


0.770 





.81 


1.11 + 0.35 


0.99 ± 0.22 


0.15 


14.43 


1803+784 


0.679 





.88 


2.79 ± 0.30 


0.99 + 0.22 


0.26 + 0.03 


13.43 


1807+698 


0.051 





.87 


1.71 + 0.32 


7.85 ± 2.44 


0.32 


14.26 


1823+568 


0.664 





.85 


1.45 + 0.21 


0.17 


0.42 


13.65 


2005-489 


0.071 





.71 


1.21 + 0.02 


9.85 + 1.71 


4.12 ± 1.77 


14.86 


2007+777 


0.342 





.91 


1.72 + 0.41 


1.17 + 0.18 


0.17 


13.66 


2131-021 


0.557? 





.96 


1.84 + 0.31 


0.16 + 0.04 


0.05 


13.16 


2200+420 


0.069 





.85 


2.14 + 0.66 


8.65 + 4.62 


0.88 


14.25 


2240-260 


0.774 





.89 


1.03 


0.26 + 0.10 


0.07 


13.32 


2254+074 


0.190 





.88 


0.56 + 0.27 


0.60 + 0.19 


0.09 


13.25 



* These data are taken from Tables 1 and 2 of Sambruna et al.,(1996). 
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Table 3. Linear Regression Analysis Results (y = Ax ± B,N = 52) 



7 



y 


X 


A 


B 


r 


P 




log L r 


O^rx 


6.64 ± 0.49 


38.24 ±0.37 


0.886 


2.26 x 10" 


17 


logL r 


Vp 


-0.67 ±0.06 


53.13 ±0.92 


-0.841 


4.42 x 10" 


14 


logic 


Oi rx 


2.45 ± 0.45 


43.94 ± 0.33 


0.614 


2.70 x 10" 


-6 


logL 


Vp 


-0.28 ±0.05 


49.93 ± 0.68 


-0.660 


2.52 x 10" 


-7 


log Lj; 


Ctrx 


-1.41 ±0.43 


46.08 ±0.32 


-0.42 


0.002 




log L x 


V p 


0.11 ±0.05 


43.43 ± 0.73 


0.30 


0.03 





r denotes the correlation coefficient and P denotes the chance probability. 
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Fig. 1. a: Plot of logL r — a rx for RBL and XBL samples. The empty circles represent RBLs, solid squares represent 
1652+398 and 2005-489, and the solid circles represent XBLs. The horizontal solid line is a rx = 0.75 and the vertical 
solid line is logL r = 43.25. The oblique solid line is the regression line logL r = (6.64 ± 0A9)a rx + (38.24 ± 0.37); b: 
Plot of log L r — log v p for RBL and XBL samples. The empty circles represent RBLs, solid squares represent 1652+398 
and 2005-489, and the solid circles represent XBLs. The horizontal solid line is \ogv p — 14.7 and the vertical solid line 
is \ogL r = 43.25. The oblique solid line is the regression line logL,. = -(0.67 ± 0.06) \ogv p + (53.13 ± 0.92). 
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Fig. 2. a: Plot of logL — a IX for RBL and XBL samples. The empty circles represent RBLs , solid squares represent 
1652+398 and 2005-489 and the solid circles represent XBLs. The horizontal solid line is a rx = 0.75; b: plot of 
logL — logi'p for RBL and XBL samples. The empty circles represent RBLs, solid squares represent 1652+398 and 
2005-489, and the solid circles represent XBLs. The horizontal solid line is logv p — 14.7. 
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Fig. 3. a: plot of logL x — a IX for RBL and XBL samples. The empty circles represent RBLs, solid squares represent 
1652+398 and 2005-489 and the solid circles represent XBLs. The horizontal solid line is a rx = 0.75; b: plot of 
logL x — \ogVp for RBL and XBL samples. The empty circles represent RBLs, solid squares represent 1652+398 and 
2005-489, and the solid circles represent XBLs. The horizontal solid line is logv p — 14.7. 



